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Abstract 

Two potential sources of tuning exist in composite Higgs models: one comes from 
keeping the Higgs VEV below the compositeness scale and one comes from keeping 
the Higgs light after constraints on the top partner masses are applied. We construct a 
measure that determines whether these tunings are independent or not and combines 
them appropriately. We perform a comprehensive scan of the parameter space for 
three explicit models and report the minimum tuning values compatible with existing 
collider constraints. Tuning values are given as functions of resonance masses and 
deviations to the Higgs couplings so the effect of future constraints can easily be 
quantified. The current minimum tuning in the minimal model is 2.5-5% and will be 
decreased to around 0.8-3.3% if no top partners are observed over the lifetime of the 
LHC. 
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1 Introduction 


For a long time the concept of naturalness has provided strong motivation to expect new 
physics at the TeV scale. As yet, evidence for this new physics stubbornly remains unforth¬ 
coming and many extensions to the Standard Model are coming under increasing pressure. 
Whether they remain natural or not is becoming an ever more common question. 

A particularly compelling way of extending the Standard Model is to replace the ele¬ 
mentary Higgs with a composite one, specifically a pseudo Nambu-Goldstone boson of a 
new, spontaneously broken, global symmetry |l]-3|. Not only is this an intriguing prospect 
by itself, given that we already know of several examples of compositeness in nature, it is 
also able to explain the fermion mass hierarchies [4j[5j and, of course, render the Standard 
Model natural, many composite Higgs sectors being insensitive to physics at scales above 
the compositeness scale. Regrettably the dearth of hints for new physics points towards 
a compositeness scale significantly higher than the electroweak scale, thus protection from 
physics at higher scales is not worth what it once was and most models reproducing the 
correct electroweak scale are expected to be tuned |6|j^] The main aims of this paper are to 
find out just how tuned composite Higgs models are now, and how this will change after 
future collider experiments. 

Several factors contribute to the tuning in composite Higgs models. The most obvious 
is the tuning required to keep the Higgs VEV below the compositeness scale. Even then it 
has been known for a while now that, in the most commonly considered models, the Higgs 
mass is further correlated with the mass of the lightest top partner |8 -14 , top partners 
being the coloured, composite fermions that allow the top quark to couple strongly to the 
Higgs. A consequence key to our study is that collider constraints on top partner masses 
cause a second tuning in the Higgs mass. To properly quantify the overall tuning in a 
given model we must therefore devise a tuning measure that determines whether any two 
tunings are independent, then combines them in a way that accounts for this. 

The other major obstacle in answering our central questions comes from the number of 
parameters. In the simplest, viable model there are still nine parameters and, by design, 
we are searching for tuned regions of the parameter space. The machinations taking a 
parameter point to a Higgs VEV and mass are not easily invertible so we could be left with 
a costly scanning problem, especially if we want to be convinced that we have sufficient 
coverage of the parameter space for the tuning values we find to be representative. Our 
strategy here is to use the nested sampling algorithm, as implemented in the Multinest 
software package 15-17 . This technique is particularly well suited to the problem but, to 


Muned composite Higgs models may still be of interest for other phenomenological reasons [7 . 
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the best of our knowledge, has not been applied to composite Higgs models before now. 

Our results comprise a set of overall tuning values as functions of phenomenologically 
interesting parameters: top partner masses, charged vector-boson resonance masses, devia¬ 
tions to the Higgs couplings, and the compositeness scale. We derive these for three explicit 
models. All are based on the minimal viable spontaneous symmetry breaking pattern in 
the composite Higgs sector, 50(5) —> 50(4) fl8|, and on the simple, 4D constructions 


developed in refs. 19,20 . We consider three different choices for the 50(5) embeddings of 
the top quark. The 5-5 model is the simplest, the 14-14 model is a slightly more involved 
model but allows for lower tuning, and the 14-1 model accommodates a fully-composite 
right-handed top quark. 

From this data we find minimal values for the tuning compatible with current collider 
constraints on the top partner masses and deviations to the Higgs couplings, and we project 
how the tuning will worsen if no new physics is seen in the imminent future. We find values 
of 3.7%, 5% and 2.5% for the current tuning in the 5-5, 14-14 and 14-1 models. Assuming 
a reach no greater than 2 TeV for top partner exclusion after 300 fb _1 at the LHC 21,22 
these values decrease to around 1%, 3.3% and 0.8% respectively. To beat the limits coming 


from top partner searches the Higgs couplings would need to be measured to an accuracy 
of a few percent. Charged vector-boson resonance searches have a much weaker effect on 
the tuning. In all three models the double tuning associated with simultaneously getting 
a light enough Higgs and the correct VEV is the main contribution. 

The rest of this paper is organised as follows. In section[2]we provide a general discussion 
on tuning in models predicting values for more than one observable, deriving a new tuning 
measure to properly account for multiple tunings. Section [3] contains an overview of the 
composite Higgs models we will be studying. Section [4] details our scanning procedure 
and our results are presented in section [5j We conclude in section [6j The appendices give 
further information on the models studied. 


2 Tuning in more than one observable 

In any model predicting a value for more than one observable there is the possibility of 
multiple, independent occurrences of tuning. For the models we will study here values are 
predicted for three relevant observables: the Higgs VEV, the Higgs mass and the mass 
of the top quark. A tuning can be defined for each and the overall tuning in the model 
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should be a combination of all threeJ3 To determine how to combine them one should first 
determine whether they are independent or not. 

If all observed tuning stems from the Higgs VEV alone being highly sensitive to the 
input parameters then the tunings are not independent and the model only really suffers 
from a single tuning. The other two observables may, for example, be parameterised in 
terms of the Higgs VEV as rrih = Xv and m t = y t v , so this kind of tuning occurs when the 
values of A and y t do not demonstrate a similar sensitivity. The other extreme scenario, 
triple tuning, occurs when all three values - v, X and yt - are highly sensitive to the input 
parameters. In other words one must tune the input parameters to get an acceptable value 
for the Higgs VEV, tune them again to get an acceptable value for the Higgs mass, then 
tune a third and final time to get an acceptable value for the mass of the top quark. Such 
a model is much more finely tuned and this should be reflected in the measure used to 
quantify the tuning. 


The most widely used measure of a single tuning is the Barbieri-Giudice tuning 23 


Abg = 


<91nO 


31nx 


( 2 . 1 ) 


O—Oe 


which quantifies the tuning of the observable O, taking the experimentally measured value 
O exp , with respect to the input parameter x. In models with several input parameters, x, 
it is often generalised to 

9ln ° ( 2 . 2 ) 


Abg = max 


d In x. 


O — O a 


One way of understanding this tuning measure is to think of the experimental result, 
O = O exp , as defining a surface in the parameter space. The Barbieri-Giudice tuning is 


then closely related to the magnitude of the normal of this surface 24,25 
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BG 
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d In O 


d In Xi 


= |V,0| D „ 


0 = 0 . 


(2.3) 


where [V/]j = d In /din Xi is a gradient operator in the corresponding logarithmic parameter 
space, In x. A large value results from a highly curved surface, matching our intuition about 
how a tuning measure ought to behave. Because of this we can think of the normals as 
defining ‘tuning vectors’ for the model. Tuning vectors have the added advantage that 

2 In a strict sense we do not actually predict the Higgs VEV as we rescale the compositeness scale to 
give the correct value. Given a full, UV complete theory the compositeness scale could not be freely scaled 
like this therefore any tuning in the ratio £ = v 2 / f 2 should still be taken into account. 
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their magnitudes are insensitive to the basis chosen for the parameter space so the tuning 
cannot be removed by making a clever choice of variables. 

When trying to evaluate the overall tuning in a set of no observables, O, we start by 
writing down individual, single tunings for each observable using the tuning vectors above 


A? = |V,0“| 0 . 0 „ p (2.4) 

where a runs over all observables and O = O exp means that all observables are set to their 
experimentally measured values. The total single tuning is then given by the sum of the 
individual contributions normalised by the number of observables 


Ar 


no 


E A t 


(2.5) 


In models that only suffer from tuning in a single observable the surfaces O a = 0“ xp will 
all be approximately aligned and the tuning vectors approximately parallel. For example, 
the extreme case for the models studied here would be for the Higgs VEV to be the only 
value that depends on the input parameters at all. We would then find 


Wim h = Vim t = V t v 


( 2 . 6 ) 


as dX/dxi = dyt/dxi = 0, and the overall tuning would be given by 


A x = - (A” + A™' 1 + A™*) = A*. 

O 


(2.7) 


This explains our choice of normalisation; we only want to count the unique source of 
tuning once in this limiting case. 

When two of the tunings are independent it seems reasonable that the individual tun¬ 
ings defined above should be multiplied together somehow. A quantity that captures this 
behaviour is 



V,O a .V / O a V ; O a .V ; 0 6 
WiO a .ViO b \/ l O b .\/ l O b 


( 2 . 8 ) 


For fully independent tunings the individual surfaces O a = 0“ xp will not be aligned at all, 
the tuning vectors will be orthogonal, and V/O a .V;O fc will vanish to give A^ 6 = A“Aj. For 
fully dependent tunings, on the other hand, the tuning vectors are parallel so V;0“.V/0 6 = 
A“Aj and Aif vanishes instead. To get an overall double tuning we evaluate this quantity 
for each possible pair of observables and take the sum. In the models studied here this 
gives 


A 2 


1 

2 



(v,m h ) 

2 




_|_ m h ,m t ) 


(2.9) 
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The normalisation is again chosen for the limiting, non-trivial case, where two tunings are 
independent and the third is fully dependent. Two terms in the sum will then be equal 
and the third will vanish to give A 2 = (for example) so that we only count the 

unique source of double tuning once. 

Finally, it is possible that all three tunings are independent and the model suffers from 
triple tuning. We can quantify this by extending the double tuning measure as defined in 


(2.8) in the obvious way 


A abc _ 

.3 ~ 


X7iO a .ViO a 

v*o a .v,o 6 

v*o a .v,o c 


WiO a .\/iO b ViO a .ViO c 

v,o ft .v,o c v,o c .v,o c 


1/2 


0 = Oe 


( 2 . 10 ) 


This quantity again behaves as desired: it is zero if any two of the tuning vectors are parallel 
(i.e. the associated tunings are dependent) and evaluates to A'j'A^’A'j; if all three tuning 
vectors are orthogonal (i.e. the tunings are fully independent). Only one such quantity 
exists in the models studied here as there are only three observables being considered. 
More generally, one would sum over all distinct combinations of a, b and c and normalise 
accordingly. 

Combining all three tuning measures gives the overall tuning 


A — Ai T A 2 T Ag 


( 2 . 11 ) 


When significant triple tuning is present in a model it will dominate to give A « A 3 . When 
significant double tuning is present but triple tuning is not one finds A ps A 2 . When only 
single tuning is present one finds A « Aj and the tuning is closely related to the more 
commonly used Barbieri-Giudice measure. 

Extending our measure to models predicting values for more observables is straight¬ 
forward. One defines the no x no tuning matrix [M/±] ab = V/0“.V;0 6 . The nth tuning 
measure between a particular choice of n observables is found by evaluating the deter¬ 
minant of the n x n submatrix containing only those rows and columns. The total nth 
tuning measure is the sum over all such submatrices normalised by no — (n — 1). This 
normalisation ensures that the one unique nth tuning is only counted once in the limiting 
case of exactly n independent individual tunings. 

Before moving on we point out that many properties of our tuning measure may well be 


captured by other measures in the literature 26,27 , particularly those based on a Bayesian 


approach 28-32 . Indeed, when there are exactly three observables the matrix appearing 


in our definition of triple tuning is the square of the Jacobian matrix encountered in a 
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Bayesian approach. It seems quite possible that our tuning measure can be rigorously 
derived using a Bayesian approach with a suitable choice of priors for the complete set of 
parameters and experimental constraints for the complete set of observables (i.e. including 
the composite sector resonance masses and so on). 


3 Model overview 


The most promising way for a composite Higgs to be realised is as a pseudo Nambu- 
Goldstone boson (pNGB) that is part of a hitherto unobserved composite sector, the com¬ 
posite sector itself emerging from a confining gauge theory (examples of suitable theories 
can be found in refs. 


33-421). Within this framework the minimal model compatible 


with custodial symmetry in the Higgs sector, and therefore precision electroweak measure¬ 
ments [43], assumes a spontaneous symmetry breaking pattern SO( 5) —>■ S'0(4) in the 
composite Higgs sector 18 . Four pNGBs are associated with this symmetry breaking 


pattern, exactly the right number to make up a Higgs doublet. Interactions of the pNGBs 
are determined by low energy theorems so all that is left is to couple the elementary, Stan¬ 
dard Model (SM) fermions to the composite sector. This is invariably achieved using the 
idea of partial compositeness [4,5], whereupon the forms of the couplings are fixed once it 
is decided which SO( 5) representations the composite operators coupling to the elemen¬ 
tary fermions are in or, equivalently, which representations the elementary fermions are 
embedded in. 

Throughout this work we will focus on 4D models like those presented in ref. fl9], also 


utilising specific constructions found in refs. 11,44 . Expressions for all important quanti 


ties are taken directly from these references. These models provide an effective description 
of the pNGBs and the lowest-lying vector-boson and fermion resonances emerging from the 
composite sector. They have the advantage of being simple, calculable and, most impor¬ 
tantly, having clear links with collider phenomenology. All of the models we will consider 
are based on the minimal symmetry breaking pattern SO(5) x U(l)x —> SO(4) x U(l)x 
(the additional U(l)x is included to give the correct hypercharge for the fermions). Since 
these models have been comprehensively explained many times in the literature we will 
only give a brief overview here. Further details are provided in appendix [A] 

We will study three embeddings for the top quarks in detail. Embeddings for the lighter 
fermions are not generally specified as their contributions to the Higgs potential tend to 
be subdominant and will be neglected]^] although we will consider modifications to the hbb 


3 An exception in some models, that we do not study here, is the right-handed tau 45 . 
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coupling as it is phenomenologically interesting. The three models and their SO(5) xU(l)x 
embeddings for the left and right-handed top quarks, q and t c , are 


The 5-5 model with q,t c G 5 2 /3 19 46 


The 14-14 model with q, t c G 14 2 /3 11,44 


The 14-1 model with q G 14 2 /3 and t c G l 2 /3 11,44 


Our choice of embeddings is strongly motivated by ref. 11 , where the authors show that 


these three models capture three qualitatively different scenarios for tuning in composite 
Higgs models. The 5-5 model relies on a cancellation between fermion contributions of 
different orders to generate the Higgs potential, the 14-14 model generates the potential 
using fermion contributions of a single order, and the 14-1 model works the same as the 
14-14 model but also allows for a fully-composite right-handed top quark. 

For each model we take the full set of composite sector input parameters - a collection 
of masses, couplings and mixing parameters - and use them to evaluate the Higgs VEV 
and the masses of the Higgs and top quark. In practise the Higgs VEV only appears in 
the ratio 

7> 2 

(3.1) 


£ EE — 

? P 


where / is the compositeness scale so we can simply rescale / to give the correct Higgs 
VEV instead of treating / as a separate input parameter. After doing this rescaling only 
points reproducing the observed values for the remaining two quantities are kept. 

For each viable point we then evaluate the spectrum of resonances predicted. Of partic¬ 
ular phenomenological importance are the masses and SM quantum numbers of the lightest 
top partners. These are coloured, vector-like fermions that mix with the top quark and 
enable it to couple to the composite Higgs. They are vital for generating both a large top 
quark Yukawa coupling and a viable Higgs potential [8-14 , and usually provide the best 
collider signal]^] Also important are the lightest vector-boson resonances charged under 
the electroweak gauge group, although these provide weaker and more model-dependent 
constraints than the top partners. 

Lastly we evaluate any deviations of the Higgs couplings to other SM fields relative to 
their SM values. These are parameterised using the ratios 


c(hxx) 


r x = 


csmPxx) 


(3.2) 


4 Similar states exist for the other SM fermions but, since the other SM fermions couple much less 
strongly to the Higgs, they have a much smaller effect on the Higgs potential and their masses are not as 
well constrained by such considerations. 
















for couplings c and csm in the composite Higgs model and SM respectively. 


3.1 Universal features 


Several aspects of the models we are interested in, particularly those concerning the gauge 
sector, only depend on the global symmetry breaking pattern so are equivalent for all 
choices of matter embedding. 

The gauge sector is described by an angle, tg = tan#, quantifying the amount of 
elementary-composite mixing in the gauge sector and taken to be small; a mass, m p for the 
lightest vector-boson resonance; and a mass, m a of another, heavier vector-boson resonance. 
We allow these parameters to freely vary in the intervals 

to G [0,1] m p ,m a G [0.5,10] TeV (3.3) 

with m a > m p . Here and in the fermion sector we check that the value of / found for each 
point is consistent with all dimensionful parameters having magnitudes less than 47 rf. 

The spectrum of massive vector-boson resonances coming from the composite sector 
includes several states charged under SU(2)l x U(1)y- The quantum numbers and masses 
(up to small, post clectroweak symmetry breaking corrections) of the lightest charged states 
are given by 

• l±i with mass m Pl = m p 


• 3 0 with mass m P3 = m p /cg 


where eg = cos 6. 

The modification to the hVV coupling, where V is a W or Z boson, is given by 

r v = \/W- (3-4) 


Modifications to the htt and hbb couplings depend on the embeddings chosen for the quarks 
and will be described in the following subsections. Up to small contributions from lighter 
states the modification to the loop-induced hgg coupling, r g , is the same as that of the 
htt coupling. The modification to the loop-induced /177 coupling is determined from the 
previous quantities via the formula 


_ A x r v + \A 1/2 r t 
r ^~ 17+ | A l/2 

neglecting contributions from states lighter than the W and Z bosons, 
standard W and Z boson and top quark loop functions, evaluating to A 
H 1/2 ~ 1-375 respectively. 


(3.5) 

The H’s are the 
1 ~ —8.324 and 
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3.2 5-5 model overview 


The simplest model embeds both the left and right-handed top quarks in the 52/3 represen¬ 
tation of SO( 5) x U( 1).y- The parameters of the composite sector comprise a Yukawa-like 
coupling, Y , coupling the top partners to the Higgs; top partner mass terms uiq, my and 
m Y ; and elementary-composite mixing parameters, d q and dt. We allow the parameters to 
freely vary in the intervals 


mQ,m T ,m Y G [0.5,10] TeV Y G [—10,10] TeV 


d q ,d t e[ 0,1]. (3.6) 


There are two top partners, also in the 52/3 representation. Breaking these down into 
SU(2)l x U(l)y multiplets the lightest states are 


• l 2 / 3 = r 2 / 3 with mass m l2/3 

• 2 1/ 6 = (T 2 / 3 , B_ i / 3 ) with mass m 2l/6 

• 2 7/6 = (T 5 / 3 , B 2/3 ) with mass m 2r/6 


The remaining modifications to the Higgs couplings are given by 

l- 2 £ 


r„ = r t = r b = 




(3.7) 


In this model the Higgs quartic coupling is generated at a higher order than the 
quadratic coupling (specifically at quartic rather than quadratic order in the elementary- 
composite fermion mixing parameters). A high degree of tuning is therefore expected. 


Further details for this model are given in appendix A.2 


3.3 14-14 model overview 

This model embeds both the left and right-handed top quarks in the 14 2 /3 representation 
of SO(5) x U(l)x- The parameters of the composite sector comprise two Yukawa-like 
couplings, Y\ and Y 2 , coupling the top partners to the Higgs; top partner mass terms 
rr iq, rriT and m Y ; and elementary-composite mixing parameters, d q and d t . We allow the 
parameters to freely vary in the intervals 

rriQ,mT,'m Y G [0.5,10] TeV Y b ,Y 2 G [—10,10] TeV d qi d t G [0,1]. (3.8) 

There are two top partners, also in the 14 2 /3 representation. Breaking these down into 
SU{2)l x U (1)y multiplets the lightest states are 
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• l 2 / 3 = T 2 / 3 with mass rai 2/3 

• 2 1 / 6 = (T 2 / 3 , B_ 1/3 ) with mass m 2l/e 

• 2 7/6 = (T 5/3 , . 62 / 3 ) with mass m 27/6 

• 3 5 / 3 + 3 2 /3 + 3_i / 3 with mass m 3 

The remaining modifications to the Higgs couplings are given by 

5(1 - 8 £ + 8 e 2 )Hi - 2(4 - 23£ + 20£ 2 )y 2 1 ~ 2£ 

rg n 2£(1 — £)[5(2£ — l)li + 2(4 — 5^)1^] n vT^' ( ' } 

Unlike the 5-5 model, the entire Higgs potential is generated at the same order (quadratic 
order in the elementary-composite fermion mixing parameters). Hence the tuning can be 
expected to be less severe. Further details for this model are given in appendix |A.3[ 

3.4 14-1 model overview 

This model embeds the left-handed top quark in the 14 2 / 3 representation of SO(5) x 
U(l)x and the right-handed top quark in the l 2 /3 representation. The parameters of the 
composite sector comprise two Yukawa-like couplings, Y x and Y 2 , coupling the top partners 
to the Higgs; top partner mass terms rriQ l , mg 2 and my; and elementary-composite mixing 
parameters, d q and A. We allow the parameters to freely vary in the intervals 

mg 1 ,mg 2 ,my, A G [0.5,10] TeV Y \, Y 2 G [—10,10] TeV d q G [0,1]. (3.10) 

There are two top partners, also in the 14 2 / 3 representation. The second top partner 
is needed to break an accidental symmetry that would otherwise increase the tuning. 
Breaking these down into SU{2)l x U(1)y multiplets the lightest states are the same as in 
the 14-14 model, i.e. 

• I2/3 = I2/3 with mass mi 2/3 

• 21 /6 = (T 2/ 3, B_ 1/3) with mass m 2l/6 

• 2 7/6 = (T 5 / 3 , B2/3) with mass ?n 27/6 

• 3 5 / 3 + 3 2 / 3 + 3 _i /3 with mass m 3 
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The remaining modifications to the Higgs couplings are the same as in the 5-5- model and 
are given by 


r g = r t = r b = 


l-2£ 


(3.11) 


Once again the entire Higgs potential is generated at the same order (quadratic order 
in the elementary-composite fermion mixing parameters) so low tuning could be expected. 
Embedding the right-handed top quark in the 1 of SO( 5) also allows for it to be fully 
composite, as there is no need for its couplings to explicitly break SO( 5). Further details 


for this model are given in appendix A.4 


4 Scanning procedure 


Given a set of input parameters, x, we wish to obtain the region of the parameter space 
in which the masses predicted for the Higgs and the top quark closely match the observed 
values. As the expressions involved cannot easily be algebraically inverted the simplest 
approach would be to scan over each parameter using a flat grid scan. Unfortunately this 
scales badly with the dimension of the parameter space, which is between nine and ten 
for the models studied here. On the other hand the nested sampling algorithm, as im¬ 
plemented in the Multinest software package, has proven particularly useful for sampling 
the non-trivial (and possibly multimodal) functions encountered in many particle physics 
and cosmology examples. We here briefly summarise the application of the technique to 
our particular problem leaving the finer details of nested sampling itself to the original 
papers 15 m■ 

Given O = {mh,rn t } the likelihood of any particular model with N p parameters x is 


P(0|x) = fjexp 


[O a 


— O a I 2 

^expJ 


2(er a ) 2 


(4.1) 


where O a (x) is the predicted value of the ith observable with experimentally measured 
value 0“ xp , a a is the error in 0“ xp , and the product runs over all observables. For our 
purposes cr a characterises how close we want the masses to be to their observed values. 
Given a prior knowledge, p(x), of the distribution of model parameters we can determine 
the posterior probability of x via Bayes’ theorem 


p(x|0) = P(0|X > P(X) . 

Zj 


(4.2) 
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The points giving higher posterior probabilities predict more viable top quark and Higgs 

masses. 

The normalisation constant, Z, is the Bayesian evidence 


Z = j p(0|x)(x)p(x)d 7Vp x. 


(4.3) 


Nested sampling is a Monte Carlo method that calculates the evidence by transforming 
this difficult, multi-dimensional evidence integral into a one-dimensional integral that is 
easy to evaluate numerically. As a by-product one obtains posterior samples and it is these 
sample points that we will interpret in the following sections. 

We use a flat prior on all parameters in this study. The choice of prior ultimately 
determines the sampling density in the final regions of interest but does not have a dominant 
effect on whether those regions are found or not. Given that the object of this study is to 
locate viable regions of parameter space, irrespective of the final sampling density, a study 
with one prior is sufficient. We have checked that doubling the number of live points in 
the Multinest algorithm from 4000 to 8000 gives similar results in terms of the located 
regions, suggesting that we have good coverage of the parameter space. 

For the Higgs and top quark masses we choose central values of 125 GeV and 155 GeV, 
and ‘errors’ of 5 GeV and 15 GeV. These choices result in points with G [120,130] 
GeV and rn t G [140,170] GeV being favoured. The low value and large interval for the top 
quark mass is chosen to account for running down to the electroweak scale and converting 
to the pole mass, the original prediction giving the running mass at a variable scale / ~ 
few TeV. Recall that the Higgs VEV only appears in the ratio £ = v 2 /f 2 so no additional 
input is required; we can simply rescale / to give the correct Higgs VEV. 

To account for collider constraints we assign a likelihood of zero to any point where 
at least one of the top partner masses does not satisfy the CMS limits m 2l/6 > 786 GeV, 
mi 2/3 > 696 GeV, m 2?/fi > 800 GeV and ms > 800 GeV 


‘ 7/6 


47 


Note that these limits only apply to a single top partner in iso’ 
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(see also refs. 


50 


ation. Since the mode 


52 


s we 


consider all contain more than one top partner the actual limits will be stronger. A full 
reinterpretation of the analyses would be required to quantify this and is beyond the scope 
of this work. Limits on the vector-boson resonance masses are not applied. In a simplified 


approach the current limits are around 1.5 TeV 53 but they are often weakened by model 


dependent effects. We anyway find that these masses are not so strongly correlated with 
the tuning so such limits are unlikely to have a significant impact on our results. Limits on 
these models from Higgs coupling measurements have been derived in ref. [54] and imply 
that / > 700 GeV or, equivalently, £ < 0.12. Constraints from electroweak precision 
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measurements and flavour-changing effects are all model dependent and are not applied 
here. Our results therefore provide a conservative lower bound on the tuning]^] 


5 Results 

After performing the scan detailed above we evaluated the tuning for each point as detailed 
in section [2j Tunings for the three models are then plotted against the mass, m p , of the 
lightest vector-boson resonance, which has SM quantum numbers l±i; the mass, my, of 
the lightest top partner resonance, which has different SM quantum numbers detailed by 
colour coding of the plots; the ratios, r x , of the Higgs couplings to their SM values, and 
the ratio, £ = v 2 /f 2 , of the Higgs VEV to the compositeness scale. These plots are shown 


in figures 5.1 5.2 and A3 In all plots we provide a simple lower bound on the tuning by 
connecting the extremal points using a convex hull. 

Double tuning dominates in all models, specifically a double tuning between the Higgs 
VEV and mass. This can be seen in the plots of tuning against £ where the lower bound 
from the single tuning, Ai, is shown alongside the overall tuning (we also include the 
naive tuning measure, l/£, in these plots for comparison). We can understand this from 
the arguments given in refs. [8-13]. A light Higgs can only be obtained without double 


tuning in the presence of light top partners but the CMS limits are already disfavouring 
this possibility. Triple tuning was not significant in any of the models we studied. 


5.1 5-5 model results 


Full results for the 5-5 model are shown in figure 5.1 The minimum tuning we find is 


A = 27, i.e. about 3.7%. Previous estimates of the tuning in this model have been a little 
less severe as they did not consider the double tuning required to get a light enough Higgs. 


This is shown explicitly in figure 5T where the single tuning is shown to be consistently 
below the overall tuning by a factor of two to three. The naive, l/£ estimate of tuning is 
a little smaller still. 

Points with milder tuning predict lighter top partners in the 27/6 SM multiplct, a lower 
compositeness scale, and a high degree of elementary-composite mixing in the top quark 


5 Additional couplings can be included that only really influence electroweak precision observables and 
the flavour sector. The extra constraints will generally impose limits, and perhaps tunings, on this extended 
parameter space while leaving the minimal parameter space largely unaffected. A caveat to this statement 
comes from the lepton sector. Increased levels of tau compositeness, motivated by flavour considerations, 


can actually lessen the dependence of the Higgs mass on the top partner masses and reduce the tuning 45 . 
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Figure 5.1: Tuning in the 5-5 model as a function of p mass, top partner mass, Higgs 


couplings and £. Further details are given in the text. 
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sector. This all leads to larger modifications to the couplings of the Higgs to gluons and 
fermions, and the tuning can be well constrained in this model by measuring these couplings 
more precisely. Once a precision better than 5% is achieved the model very quickly moves 
beyond A = 100 and begins to look unnatural. Further constraints on the top partner 
masses also quickly lead to a higher degree of tuning. With a final top partner reach of 
about 2 TeV the LHC can push the tuning to A > 100 from these searches alone. 

Points giving the correct Higgs mass in this model clearly separate into two regions: 
one in which the 27/6 SM multiplet tends to be the lightest top partner and one in which 
it tends to be the I 2 / 3 . The latter region is significantly more tuned. This is because it is 
harder to keep the I 2/3 light than the 2 7 / 6 . The I 2/3 mixes directly with the elementary, 
right-handed top quark so its mass is more constrained by the observed top quark mass. 
The 2 7/0 does not mix directly with the elementary top quark so its mass is more flexible. 


5.2 14-14 model results 


Full results for the 14-14 model are shown in figure 5.2 The tuning is a little less severe 
than in the 5-5 model, particularly for higher top partner masses, and the minimum we find 
is A = 20, i.e. 5%. The 14-14 models tend to have a much higher compositeness scale and, 
therefore, much more SM-like Higgs couplings (note that 7 y % r t in this model and only the 
minimum tuning for is shown). Lighter top partners in the 2 7 / 6 SM multiplet are still 
preferred suggesting that the composite sector is not too strongly coupled: m T 9rf so 
the composite sector coupling, gx, cannot be too large. A very high degree of elementary- 
composite mixing in the top quark sector is preferred, at the limit of the range of viability 
in our effective theory approach. Interestingly the naive, l/£ estimate of tuning actually 
overestimates the tuning here. 

The tuning can be most constrained in this model by improving limits on the top 
partner mass. Even then the change in the constraint is mild; a limit A > 30 is reached 
for top partners heavier than 2 TeV. In the Higgs coupling sector a precision better than 
2% is required to outperform this. Hence the 14-14 model will remain relatively natural 
for the foreseeable future. 


5.3 14-1 model results 


Full results for the 14-1 model are shown in figure 5.3 The tuning behaves similarly to 
that of the 5-5 model, albeit starting from a slightly higher base. The minimum we find is 
A = 40, i.e. 2.5%. The compositeness scale tends to be lower than in the 14-14 model so 
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Figure 5.2: Tuning in the 14-14 model as a function of p mass, top partner mass, Higgs 
couplings and £. Further details are given in the text. 
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constraints from Higgs coupling measurements are more important. Even so the best way 
to constrain the tuning in this model is via top partner searches, specifically for those in 
the 2 7 / 6 SM multiplet, and excluding top partners lighter than 2 TeV pushes the tuning 
to A > 130. A precision better than 3% in Higgs coupling measurements is required to 
outperform this. 


5.4 Discussion 


A common feature of all models studied here is that the severity of the tuning is being 
driven by constraints on the top partner masses, for the immediate future at least. After 
300 fb _1 of 14 TeV collisions at the LHC top partners lighter than 2 TeV are expected 
to be excluded 21,22 . Comparing this with the expected precision to which the Higgs 
couplings can be measured with the same data set, around 9% for hWW and worse for the 


other couplings 55 , and it is clear that searching for top partners is a much more powerful 


way to probe naturalness. 

In the more distant future this may change. A higher energy proton-proton collider 
could constrain top partner masses even further; lower limits of around 5 TeV and 9.5 TeV 
can be expected after 3000 fb -1 of 33 TeV and 100 TeV collisions respectively 


21 


22 


leading to tunings worse than 0.1% in the models studied here. However, similar constraints 
on the tuning can be expected from a electron-positron collider sat on the Higgs resonance. 


TLEP 56 , for example, would measure the hZZ coupling to an accuracy of about 0.15%, 


implying a tuning worse than 0.1% in the 5-5 and 14-1 models (due to its high compositeness 
scale top partner searches remain more powerful in the 14-14 model). 

The main advantage of using an electron-positron collider to probe naturalness in this 
context is that it would be able to place meaningful constraints on a wider variety of 
models. Constraints coming from top partner searches rely on a connection between the 
Higgs mass and the masses of the top partners. This holds for the models studied here, 
where the large top quark Yukawa is explained by top quark compositeness, but other 


classes of model do exist. A composite right-handed tan can change this conclusion 45 


as can the use of colourless states to stabilise the composite Higgs mass, like in composite 


twin Higgs models 57 , for example. A proper quantification of the tuning in such models 


is beyond the scope of this work, but we expect that the relationship between modifications 
to the Higgs couplings and the tuning will be qualitatively similar to the models that we 
have studied. 
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6 Conclusions 


Double tuning - i.e. independent tunings to simultaneously get the right Higgs VEV and 
mass - is important in composite Higgs models and should be accounted for. This can be 
done by using a tuning measure like the one we construct in section [2j Despite this extra 
source of tuning minimal 4D models based on an SO(5) —» SO (4) symmetry breaking 
pattern remain fairly natural, the current minimum tunings being 3.7%, 5% and 2.5% 
in the 5-5, 14-14 and 14-1 models. If top partners lighter than 2 TeV are excluded the 
tuning worsens to around 1%, 3.3% and 0.8% respectively. We arrived at these values 
after performing a comprehensive scan of the full parameter space using a nested sampling 
algorithm, then applying the latest collider constraints on top partner masses and Higgs 
coupling deviations. Other, more model dependent constraints can be applied so our values 
provide a conservative minimum. 
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A Model details 


The models we consider consist of two sites - an elementary site and a composite site. 


All SM degrees of freedom other than the Higgs live on the elementary site 19,20 . The 


composite site comprises a non-linear sigma model describing the pNGB Higgs and a set of 
vector-boson and fermion states representing the resonances. The two sites are connected 
via a set of link fields described by a second non-linear sigma model. The link fields are 
responsible for mixing between the two sites and for giving masses to the vector bosons on 
the composite site. 

Initially the elementary site has a global symmetry G e = SO (5) x 1/(1)x and the 
composite site has a gauge symmetry G c = SO( 5) x U(l)x- The gauge fields associated 
with G c are our vector-boson states, p and px- The overall G e x G c symmetry is then 
spontaneously broken to the global, diagonal G = SO( 5) x U(l)x subgroup. The NGBs 
associated with this symmetry breaking are our link fields, Q. Since G c is gauged the 
link fields are eaten by the G c gauge fields to produce massive, vector-boson resonances 
transforming in the adjoint representation of G. At the same time G c is spontaneously 
broken to an H c = SO( 4) x U(l)x subgroup on the composite site. The NGBs associated 
with the G c —y H c symmetry breaking, E, will go on to provide our Higgs. 

Upon including the SM degrees of freedom on the elementary site G is explicitly broken 
to the usual, SU(2) L x U( l)y electroweak symmetry of the standard model. This is pre¬ 
cipitated by gauging only the SU(2)l x U(1)y subgroup of G e and by embedding the SM 
fermions in incomplete representations of G e . Nonetheless, the explicit breaking is weak 
(couplings and masses on the composite site are much greater than their elementary site 
counterparts) so the symmetry structure discussed above remains approximately intact 
and the Higgs, embedded in E, remains light due to its pNGB nature. 

The Higgs is embedded into the spurion 




0 0 0 


where n = ir a X a contains the four NGBs for the four broken SO(5) generators, {X a }, 
Sh = sin (h/f) and Ch = cos (h/f). The NGB decay constant is denoted by / and we work 
in the unitary gauge such that [7T a ] = (0, 0,0, h). The spurion transforms as a 5 of SO( 5) 
and acts a convenient book-keeping device when determining the interactions of the pNGB 
Higgs allowed by the symmetries of the model. 

The Lagrangian for these models splits into contributions from the elementary site, the 
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composite site and the mixing terms. The contribution from the elementary site is 


= -4k w ^ ' + Gsb) «*•» + (sb) w + - (A - 2) 

where PF and 5 are the usual SU(2) L and U( l)y held strength tensors; q is a (two- 
component) fermion for the third-generation, left-handed quark doublet; t c is a fermion for 
the right-handed top quark, and D e is the covariant derivative involving the elementary 
gauge fields. Dots denote terms involving the lighter fermions and, for later convenience, 
the quark kinetic terms are not canonically normalised - the normalisation factor will be 
explained shortly. 

The contribution from the composite site is 

Cc = - - l^Px^Px + + Q0cQ + Q^IPcQ 0 (A.3) 

+ TiJj) c T + T c tJp c T c - m Q QQ c - m T TT c - m y QT c - Y(E)QT C + h.c. + ... 


where p and px are the held strength tensors for the composite, G c gauge held; E is the 
real scalar (A.l) containing the Higgs helds, and D c is the covariant derivative involving 
the composite gauge helds. Q and T are (two-component) fermions, each coming in a 
Dirac pair. These fermions have diagonal mass terms, mg and rn y , an off-diagonal mass 
term, my, (assuming it is consistent with G c , which it is in the models studied here) and 
Yukawa-like terms, Y(E), coupling the fermions to the Higgs via the spurion. The forms 
of the Yukawa-like terms depend on which representations of G c the fermions belong to. 
Note that Q C T terms are not present despite being allowed by all symmetries of the model. 
These terms are omitted to ensure that the Higgs potential remains finite. 

Finally, the mixing terms are 


£ m = ^f(D e+Ct ^)(D^ +c n) t + A [R q (tt)qQ c + R t {Q)t c T] + h.c. + ... (A.4) 

where Q is the complex scalar parameterising the G,, x G c G NGBs. These transform 
under both G e and G c so the covariant derivative, D e+C , contains both elementary and 
composite gauge helds. The subsequent terms then use D to mix q and t c with Q c and T 
in a way respecting the initial G e x G c symmetry (or at least they would if q and t c came 
in complete G e representations). This is done by using projections, i?(D), corresponding 
to the representations of G e that q and t c are embedded in. Owing to the non-canonical 
normalisation of the elementary fermions the actual couplings associated with the mixing 
terms for q and t c go like d q mQ and d t rriT respectively. The common scale, A, is cancelled 
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in most expressions for physical quantities, although persists as an independent variable 
in the 14-1 model in which d t is also absorbed into A 

In the end we are mostly interested in the low-energy effective theory coming from the 
above Lagrangian, which can be derived by integrating out degrees of freedom from the 
composite site: p, Q and T. Recall that the link fields are eaten by the p’s so will not appear 
in the effective theory either. The detailed structure on the composite site gets encoded 
in momentum-dependent (i.e. non-local) terms in the effective theory’s Lagrangian. In 
momentum space this Lagrangian is 

£ eff = \pI v [rVC p 2 , h)w,w u + n B (p 2 , h)B,B u + iWp 2 , h)w*B v ] + 

nq(p 2 , h)qfq + IR(p 2 , h)t c ft c + M(p 2 , h)tt c + h.c. (A.5) 

where II and M denote the form factors encoding the effects of the composite sector and 
P T is the transverse projection operator. Explicit expressions for the form factors can be 
derived once the embeddings for the elementary fermions have been chosen. 

The one-loop Higgs potential is found to be 

r°° dr ? 2 

-2 N c J J^p 2 ln [p 2 (l + U q (p 2 ,h))(l + U t (p 2 ,h)) + \M{p 2 : h)\ 2 ] (A.6) 

keeping only the contributions from the SU(2)l gauge fields and the top quark. It is more 
usually expanded in powers of Sh to give 

v(h) = -74 + Psl (A.7) 


The Higgs VEV and mass found from this potential are 


^ 2/3 


2 b/3 

m h = - p ^-0 


and (A.5) yields a top quark mass 


m t = 


M{ 0, v) 


v/n g (o,v)n t (o,i;) 


(A. 8 ) 


(A.9) 


A.l Gauge sector variables and form factors 


Form factors in the gauge sector depend only on the symmetry breaking pattern so are the 
same in all models studied here. We vary g p , f c and /q via the mixing angle and masses 


9 2,e 

9 P 


m, 


_ l r 

2 9pJ( 


™ a = 2 ^p(/c + fi 


(A. 10 ) 
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The form factor for the W boson is 


pV-(i + *2K) . 1 2 

11 g\{\ + tg)(p 2 — m 2 ) 


2 rrip{m[ 


~ ™ 2 p )t 


9%(l + t 2 e)(p 2 -™ 2 a )(p 2 ~ 


(A.ll) 


where g 2 is the observed SU(2)i gauge coupling. Plugging into (A.6) and performing the 


integral results in a contribution to the s 2 h coefficient in (A.7) 

' 2 - m 2 )t 2 


7 9 = 


9m p (m P 


647r 2 (m 2 — (1 + t 2 e )m?) 


In 


m„ 


(1 + tl)m 2 f 


(A.12) 


at leading order in tg. There is no equivalent contribution to the s\ coefficient and higher 
orders terms in Sh or tg will not be considered. 


A.2 5-5 model variables, form factors and masses 


This model was developed in detail in ref. 19 . The form factors can conveniently be 
expressed in terms of the functions 


n 5 (mi,m 2 ,m 3 ) = 


M 5 (mi,m 2 ,m 3 ) = 


(ml + nig — p 2 ) A 2 


p 4 — (m\ + m 2 + m|)p 2 + m\ml 
m\m 2 m 3 A 2 

p 4 — (m\ + ml + m|)p 2 + rrifm'l 


(A.13) 


as 


A 2 - 1 ^ 

n ? = J 2 + n s( m Q’ m T , ™-y) + 2 S h n 5 (m Q , m T , m Y + F) - U 5 (m Q , m T , m Y ) 

A 2 - r- 

n t = + n 5 (m T , m Q , m Y + F) + s 2 h II 5 (m T , m Q , m Y ) - n 5 (m T , m Q , m y + F) 


M = ShCh 


M 5 (m Q , m T , m Y + F) - M 5 (m Q , m r , m y ) 


(A.14) 


Plugging into (A.6) and performing the integral results in contributions to the and s \ 
coefficients in (A.7) at, respectively, leading (quadratic) and subleading (quartic) order in 
dt and d q . Higher orders terms in Sh or d q and d t will not be considered. The top partner 
masses, up to small, electroweak corrections, are found as follows 


• p 2 = is a zero of at Sh = 0 

• p 2 = m| i/g is a zero of I\ q at sy — 0 

• P 2 = m l 7/6 is a pole of U 5 (m Q , m T , m Y ). 
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A.3 14-14 model variables, form factors and masses 


This model was developed in detail in ref. [XT] (see also ref. |44|). The form factors can 
conveniently be expressed in terms of the functions 

(m 2 + m\ — p 2 ) A 2 


as 


fli 4 (mi, m 2 , m 3 , m 4 ) = 
Mi 4 (mi,m 2 ,m 3 ) = 


p 4 — (ml + m 2 + m 2 )p 2 + m§m| 
mpm-iVd-z^ 2 


p 4 — (mf + m 2 + m\)p 2 + m 2 m 2 


A 2 


n g = -^- + IIi 4 (mr, my + ^r, mg, m T ) + 




ni4(my, my, m Q , m T ) - 2fli 4 (mr, my + m Q , m T ) 


A Sh 


+ n 14 (m T , my + f,mQ, m T ) 
3IIi 4 (m r , my, m Q , m T ) - 8n 14 (m T , my + m Q , m T ) 


+ 5n 14 (m T , m y + , m Ql m T ) 


(A.15) 


(A.16) 


A 2 

: t = + n 14 (m Q ,my + y,mQ, m T ) + 


lli — ~^r 
a 

5 


2 ' S/l 
1,4 
16 h 


ni 4 (m Q , my + |,mg, m T ) - n 44 (m Q , m y + ^,mg, m T 
3ni 4 (mg, m y , mQ, m T ) - 8n 14 (mg, my + 1 m Q , m T ) 

+ 5fl 14 (mQ, my + ^,mQ, m T ) 


+ 


(A.17) 


,, ^ 

M = —s h c h 


2 

V5 


~ S h C h 


Mu(m Q , m T , my + ^) - M 14 (mQ, m T , my + 1 
3M 14 (mQ, m T , my) - 8M M (m Q , m T , my + ^) 


+ 5Mi 4 (mQ, m T , m y + ^) 


(A.18) 


Plugging into (A.6) and performing the integral results in contributions to both the s 2 h and 
coefficients in (A.7) at leading (quadratic) and subleading (quartic) order in d t and d q . 


Higher orders terms in Sh or d q and d t will not be considered. The top partner masses, up 
to small, electroweak corrections, are found as follows 
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• p 2 = mf is a zero of 11* at Sh = 0 

• p 2 = m | is a zero of n g at = 0 

• p 2 = m| 7/g is a pole of ni 4 (m T , my + mg, m T ) 

• p 2 = m\ is a pole of my, mg, my). 


A.4 14-1 model variables, form factors and masses 


This model was developed in detail in ref. |TT] (see also ref. [44]). There are two 14’s, Q\ 
and Qii with separate mass terms and no T as the fully-composite t c is already able to 
play this role j^] Hence dt is absorbed into A and A is treated as the independent variable. 
The form factors can conveniently be expressed in terms of the functions 


n 5 (mi,m 2 ,m 3 ) 

ni(mi,m 2 ,m 3 ) 

Mi(mi,m 2 ,m 3 ) 


(t7t| + m\ — p 2 ) A 2 
p 4 — (m 2 + m 2 + m§)p 2 + rn 2 rri 2 
(m| — p 2 ) A 2 

p 4 — (m 2 + ml + m|)p 2 + m\m 2 
77i i (m| — p 2 )A 2 

p 4 — (mf + m 2 + m 2 )p 2 + rrifrn 2 


(A.19) 


as 


A 2 - 

n 9 = + n 5 (mg 1 ,mg 2 ,my + ^) + 


d 2 q 


-Su 


— Sh 


n 5 (m Ql , mQ 2 , m Y ) - 2U 5 {m Ql ,mQ 2 ,m Y + y) 
+ n 5 (mg 11 mg a ,my+ 

3n 5 K,mg 2 ,my) - 8n 5 (777Q 1 ,777Q 2 ,my + 

+ 5n 5 (mg 1 ,mg 2 ,my + 


and 


(A.20) 


n f = 1 + n* (mg,, mg 2 , A) 

M — Mi (m Ql , m Q2 , A). (A.21) 

6 Without the second 14 an accidental symmetry causes a part of the form factor crucial for reducing 
the tuning in this model to vanish. 
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Plugging into (A.6) and performing the integral results in contributions to both the sj, and 
s\ coefficients in (A.7) at leading (quadratic) and subleading (quartic) order in d q . Higher 
orders terms in Sh or d q will not be considered. The top partner masses, up to small, 
electroweak corrections, are found as follows 


• p 2 = is a zero of n* at Sh = 0 

• p 2 = is a zero of n g at Sh = 0 

• p 2 = ml r/e is a pole of ft 5 (m Ql ,m Q2 ,m Y + y) 

• p 2 = m\ is a pole of n 5 (m Ql ,m Q2 ,m Y ). 
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